Heavy flavor suppression: role of hadronic matter 
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The role of hadronic matter in the suppression of open heavy flavored meson has been studied. 
The heavy-quarks (HQs) suppression factors have been calculated for nuclear collisions at Relativis- 
tic Heavy Ion Collider (RHIC) and Large Hadronic Collider (LHC) energies to demonstrate the 
relevance of the hadronic phase. It is found that the suppression in the hadronic phase at RHIC 
energy is around 20% — 25% whereas at the LHC it is around 10% — 12% for the D meson. In case 
of B meson the hadronic suppression is around 10% — 12% and 5% — 6% at RHIC and LHC energies 
respectively. Present study suggests that the suppression of heavy flavor in the hadronic phase is 
significant at RHIC. However, the effect of hadronic suppression at LHC is marginal, this makes the 
characterization of QGP at LHC less complicated. 
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I. INTRODUCTION 

One of the primary aims of the ongoing heavy-ion col- 
lision experiments at RHIC and LHC energies is to create 
a new phase where the bulk properties of the matter are 
governed by the light quarks and gluons, such a phase 
of matter is called Quark Gluon Plasma (QGP). In this 
context, heavy flavors (HFs) play a vital role since they 
do not constitute the bulk part of the system as they are 
Boltzmann suppressed due to their higher masses. HFs 
are produced in the early stage of the collision as their 
production is associated with large momentum transfer. 
Since HFs are not frequently created or annihilated in the 
QGP, they are able to witness the entire space-time evo- 
lution of heavy ion collision and act as a effective probe 
in QGP diagnostics. In particular the depletion of high 
transverse momentum (j>t) hadrons (D and B) produced 
in Nucleus + Nucleus collisions relative to those produced 
in proton + proton (p+p) collisions has been consid- 
ered as an indicator of QGP formation. The STAR [J, 
PHENIX @, i] and the ALICE 3 collaborations have 
measured this high px depletion. 

Several attempts have been made to diagnose the QGP 
phase using HFs [o1.l7l4l3.l26j as probes. In most of these 
studies, the role of the hadronic matter is ignored. But 
once the temperature of the system (QGP) produced at 
RHIC and LHC energies goes below the transition tem- 
perature (T c ) due to expansion, the system reverts to 
the hadronic phase, which indicate that the existence of 
the hadronic phase is inevitable. Therefore, to make the 
characterization of QGP reliable the role of the hadronic 
matter should be taken into consideration and its contri- 
bution must be disentangled from the experimental ob- 
servables. In this work an attempt has been made to 
estimate the effect of the hadronic phase on the nuclear 
suppressions of HFs. 

We studied the evolution of the HFs in the follow- 
ing scenario. We assume that the light quarks and glu- 
ons form a thermalizcd matter and the non-equilibrated 



heavy quarks (HQs) are moving through the expanding 
QGP background. While the evolution of the expand- 
ing QGP is described by the relativistic hydrodynamics 
with initial temperature and the thermalization time con- 
strained by the measured charged particle multiplicity, 
the motion of the non-equilibrated HQs arc described by 
the Fokker-Planck equation (FP) with drag and diffusion 
co-efficients arising due to the interaction of HQs with 
the expanding QGP background. The initial conditions 
for the distributions of HQs have been taken from the 
NLO pQCD results obtained for pp collisions by using 
the MNR code fijl. 



The expanding QGP converts to hadronic system when 
it cools down to the transition temperature, T c . The solu- 
tion of the FP equation for the charm and bottom quarks 
at the transition point is folded with the Peterson frag- 
mentation function [i~6j | to obtain the momentum dis- 
tributions of the heavy flavoured mesons containing the 
effects of the interaction of the expanding QGP back- 
ground. The hadronic matter evolves in space and time 
described by relativistic hydrodynamics till the matter 
gets dilute enough to freeze-out kincmatically. The mo- 
tion of the non-equilibrated HF mesons (D and B) in the 
expanding hadronic system is again described by the FP 
equations with drag and diffusion coefficients evaluated 
due to their interactions with hadronic matter composed 
of pions, kaons, eta and nucleons. The solution of the 
FP equation for the D and B mesons at the freeze-out 
point encompassing the effects of drag of both the QGP 
and the hadronic phases has been used to determine the 
suppression in the high px domain. 



In the next section a brief outline of the theoretical 
formalism used in the present study has been discussed. 
Results are presented in section III and section IV is ded- 
icated to summary and discussions. 
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II. THEORETICAL FORMALISM 

The FP equation describing the motion of the non- 
equilibrated degrees of freedo m ( dof) in the bath of the 
equilibrated dof reads as 0, US] , 



The suppression of high px D or B mesons the QGP 
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where / is the momentum distribution of the non- 
equilibrated dof, Ai and Bij are related to the drag and 
diffusion coefficients. The interaction between the probe 
and the medium enter through the drag and diffusion 
coefficients. 

We now briefly recollect the scenario for the HQs mo- 
tion through the expanding QGP. During their propa- 
gation through the QGP the HQs dissipate energy pre- 
dominantly by two processes [IHQjil]: (i) collisional, e.g. 
gQ —> gQ, qQ —t qQ and qQ — > qQ and (ii) radiative 
processes, i.e. Q + q^-Q + q + g and Q + g — >• Q + g + g ■ 
The dead cone and Landau-Pomeronchuk-Migdal (LPM) 
effects on radiative energy loss of heavy quarks have been 
considered. Both radiative and collisional processes of en- 
ergy loss are included in the effective drag and diffusion 
coefficients 
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The solutions of the FP equation have 
been convoluted with the Peterson fragmentation func- 
tions to obtain the D and B meson spectra at T c ~ 170 
McV. We omit the detailed description here as it is avail- 
able in 
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In previous works we have calculated the heavy me- 
son drag and diffusion coefficients in the hadronic matter 
consisting of pions, kaons, eta and nucleons 0, Hl[ (see 
also other recent works [22(-[2J]). The interaction of the 
D and B meson with the thermal hadrons consisting of 
pions, kaons and eta has been obtained from scattering 
lengths evaluated at NNLO usin g H eavy Meson Chiral 
Perturbation Theory {HM\PT) [21j for the elastic pro- 
cesses: D{B) + 7T -> D(B) + n, D{B) + K -> D(B) + K, 
D(B) + r) -> D(B) + 77. It is found that the D and B 
mesons interactions with the nucleons for the reactions 
D(B) + N -> D(B) + N is very small in compari- 
son to their interactions with the light mesons described 
above. 

The motion of the out-of-equilibrium HF mesons (D 
and B) in the expanding hadronic matter (HM) (com- 
posed of pions, kaons, nucleons and eta) is studied by the 
FP equations. The solution at the end of the QGP phase 
is convoluted with the Peterson fragmentation function 
to get the pt distribution of the HF mesons, which has 
been used as the input (initial ) to the FP equation for 
the hadronic phase. The drag and diffusion coefficients 
of the D and B mesons have been calculated due to their 
interactions with pions, kaons, eta and nucleons. The 
solution of the FP equation for the D and B mesons in 
the expanding HM at the freeze-out is employed to de- 
termine the nuclear suppression. The expansion of the 
background medium (either QGP or HM) is described 
by relativistic hydrodynamics 25| with equation of state 
(eos) which leads to the velocity of sound, c s — 1/ 



phase, R^a, in the QGP phase is given by: 
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where fg is given by the convolution of the solution of 
the FP equation at the end of the QGP phase with the 
HQ fragmentation to D or B meson and fi is the func- 
tion obtained from the convolution of the initial heavy 
quark momentum distribution with the HQ fragmenta- 
tion function for heavy mesons. Similarly the suppression 
factor in the hadronic phase alone can be written as, 



u aa — ~r~ 

JQ 



(3) 



fu is the solution of the FP equation describing the evo- 
lution in the hadronic phase at the freeze-out. 

The net suppression of the HFs during the entire evo- 
lution process - from the beginning of the QGP phase to 
the end of the hadronic phase is given by: 
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FIG. 1: Variation of D meson suppression at RHIC energy 
for QGP, hadronic and hadronic+QGP phase. 



III. RESULTS 

Wit hin the ambit of the model described above we have 
evaluated the suppression of the HFs separately in the 

QGP (Raa) and HM ( r aa)- Tne nct suppression is 
obtained from these two factors using Eq. [U The re- 
sults for the D meson at RHIC energy is depicted in 
Fig. [TJ We have taken the initial temperature, = 0.4 
GeV and thermalization time, n = 0.2 fm/c. These val- 
ues are constrained by the measured hadronic multiplic- 
ity, dN/dy = 1100. We observe that the D meson sup- 
pression in the hadronic phase is around 20% — 25% for 
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FIG. 2: Variation of B meson suppression at RHIC energy 
for QGP, hadronic and hadronic+QGP phase. 
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FIG. 3: Variation of D meson suppression at LHC energy for 
QGP, hadronic and hadronic+QGP phase. 



Pt = 3 — 10 GeV at RHIC energy. This suggests that 
the effects of the hadronic medium on the charmed me- 
son suppression is no n- negligible. Therefore, these effects 
should be excluded from the experimental data to esti- 
mate the suppression in QGP and make the characteriza- 
tion of QGP more definitive. Drag coefficients quantify 
the momentum degradation whereas diffusion coefficient 
dictates the momentum randomization and energy gain. 
The simultaneous effect of drag and diffusion on the dis- 
tribution shift the high px particles to the low px domain 
resulting in Raa > 1 at low pt . The results for B meson 
is displayed in Fig [5] at RHIC. In the hadronic phase the 
B meson suppression is around 10% — 12%, indicating 
greater suppression of D than B. However, the overall 
suppression of B is also less than D. Because the drag of 
b quarks (B mesons) in QGP (HM) is smaller than that 
for c quarks {D meson). 

In Fig|3]the variation of Raa against px is depicted for 
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FIG. 4: Variation of B meson suppression at LHC energy for 
QGP, hadronic and hadronic+QGP phase. 
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FIG. 5: (color on-line) The pt spectra of the D mesons at 
RHIC obtained by convoluting the charm quark to D me- 
son fragmentation with (i) the initial charm quark distribiuon 
(black) , (ii) solution of the FP equation at the transition point 
(red) and (iii) the solution of the FP equation at the end of 
the hadronic phase (i.e. at freeze-out) contains the effects of 
suppression in the QGP as well as hadronic phases (blue). 



the D meson at LHC energy. We have taken the value of 
Ti = 550 MeV and n = 0.1 fm/c for ^/JWn = 2.76 TeV. 
It is found that the D meson suppression in the hadronic 
phase at LHC energy is around 10%— 12% for pt = 3 — 10 
GeV. Fig. H] displays the depiction of B mesons at LHC in 
the hadronic as well as in the QGP phase. The effects of 
the hadronic phase is found to be negligibly small. Indi- 
cating the fact that the response of the hadronic medium 
is less pronounced at LHC than RHIC. Therefore, the 
role of hadronic medium in characterizing the QGP by 
using heavy flavours can be ignored making the task of 
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FIG. 6: Same as Fig. [5] for LHC conditions. 
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and HM. The HM is composed of pions, kaons, nucle- 
ons and eta. While the HF suppression in QGP is used 
as a signal of QGP, the hadronic suppression is treated 
as background. We observe that the suppression of D 
is more than B in the hadronic medium because the 
hadrons drag the D more than the B [2l|. The sup- 
pression at RHIC energy is non-negligible and hence the 
hadronic contributions should be taken into account in 
analyzing the experimental data. It is also interesting 
to note that the role of hadronic medium in HFs (es- 
pecially for B) suppressions at LHC is not significant 
enough. Because D and B meson distribution is harder 
at LHC than RHIC. Since the role of hadrons in B me- 
son suppression is very minimal, therefore, B may play a 
unique role in characterizing QGP. This has a great ad- 
vantage compared to other signals of QGP, f or example, 
in case of electromagnetic probes (see [26r428| for review) 
i. e. direct photons and lepton pairs the role of hadronic 
matter is significant, which makes the task of extract- 
ing QGP properties difficult after filtering out hadronic 
contributions. 



QGP detection less complex. 

The differences in the magnitude of B^ A in RHIC and 
LHC can be understood from the corresponding results 
plotted in Fig. [5] and Fig. [5] The temperature of the 
hadronic system for both RHIC and LHC varies from 
T c to T f (170 to 120 MeV) and therefore, the value of 
the drag coefficients remains same for both the cases. 
However, the input distribution to the hadronic matter 
obtained from the convolution of the HQs distribution 
at the end of the QGP phase with the corresponding 
fragmentation function is harder at LHC than RHIC , 
resulting in less suppression at LHC. 
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IV. SUMMARY AND DISCUSSIONS 

To summarize, in this work we have evaluated the sup- 
pression of HFs due to their interactions with the QGP 
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